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We study the thermodynamics of helium at densities relevant for white dwarf physics. We find 
evidence that, as the temperature is increased, there is first a first order transition between two 
superconducting phases followed by a second order transition to the normal state. These transi- 
tions occur, for realistic densities, at temperatures below the crystallization temperature and the 
crystalline state is likely to remain as the true ground state of the system. The calculations are 
performed with a screening but non-dynamical electron background and we comment on the impact 
of this and other approximations to our result. 

I. INTRODUCTION 

Helium white dwarfs (He WDs) are astrophysical objects which are composed predominantly of helium nuclei and 
degenerate electrons. At typical WD densities, the nuclei are much closer together than typical atomic sizes but 
are still widely separated compared to typical nuclear sizes. It has long been known that as WDs cool the nuclei 
crystallize, locked into position by their mutual Coulomb interactions[l-3]. Recently, it was pointed out that in He 
WDs the temperature at which the helium nuclei form a Bose-Einstein condensate (BEC) might be higher than the 
crystallization temperature and an intermediate superconducting phase may exist between the plasma and the crystal 
phases [4-9]. In this phase, it is the ions that are superconducting; the electrons form an ordinary Fermi liquid. The 
low temperature properties of this phase are dominated by the physics of an unusual "phonon" excitation [10] and 
leads to a very small specific heat and enhanced neutrino emission [11], with possible consequences for the cooling of 
He WDs [12]. A similar phase could also exist in a deuterium layer in brown dwarfs [13] and be relevant for inertial 
confinement [14-16] as well as other kinds of experiments [17, 18] where high densities are also achieved . 

It is guaranteed that at large enough densities there will be a range of temperatures where the BEC can exist while 
the Coulomb crystal cannot. This can be understood by simple scaling arguments: a BEC should form when the 
thermal de Broglie wavelength \/2n /MT (here M is the ion mass and T is the temperature) becomes comparable to 
the interparticle spacing I, so that the condensation temperature should scale as Tbec ~ 1/Ml 2 . A Coulomb crystal 
should melt when the thermal energy is comparable to the nearest-neighbor interaction, so that Tboulomb ~ Z 2 a/l, 
where Z is the atomic number of the crystallized nuclei and a = e 2 /Air ~ 1/137, where e is the size of the electron 
charge. Since the number density n ~ Z~ 3 , Tbec ~ n 2 / 3 /M while Xboulomb ~ n 1 ^ 3 . Thus, at very high density, the 
crystallization temperature is markedly lower than the condensation temperature, and for intermediate temperatures, 
the system should be a BEC. 

The natural question is: are astrophysical densities in this interesting regime? To answer this quantitative question, 
one needs to know the numerical coefficients that specify these critical temperatures. Because the condensation tem- 
perature scales inversely with the ion mass, the density at which Tbec = ^Coulomb ancl beyond which Tbec > ^Coulomb 
is smaller for lighter nuclei. Thus, if a nuclear condensate forms in WDs, it should be most easily established in He 
WDs and not carbon-oxygen WDs. 

Detailed studies have determined the crystallization temperature to be Tc ou iomb ~ (Ze) 2 / 180Z [19-21], meaning 
^Coulomb ~ (a /^)7000K, where ao is the Bohr radius. There are various suggestions for the proportionality constant 
in Tbec- Simply equating the de Broglie wavelength to the interparticle spacing suggests Tbec = 2n/Ml 2 s» 6.2/MZ 2 . 
A free Bose gas has Tbec = T 1 ] ' = 27r(47rC(3/2)/3)" 2/3 /AfP « 1.27/A/Z 2 , where £ is the Riemman zeta function. 
The temperature Tbec is expected to go up when one considers repulsive interactions [22]. A slightly more detailed 
estimate (see Ref. [6]) suggests Tbec = 4-7r 2 /3Af^ 2 13.2/MZ 2 , which is qualitatively supported by the numerical 
calculations in Ref. [23]. It is the object of this paper to make a reliable estimate of Tbec- 

A BEC composed of nuclei (and not whole atoms) with a background of degenerate electrons is a novel system with 
rich phenomenology. Because the condensed nuclei are charged, the substance is electrically superconducting. The 
electrons provide a neutralizing electric charge, and additionally the dynamical response of these electrons implies 
an unusual gapless quasiparticle[10]. These quasiparticles imbue the substance with a very small specific heat [10]. 
Moreover, these quasiparticles can annihilate into neutrinos, and the power emitted per unit volume scales like 
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T 11 , so that the phenomenological relevance of this annihilation for He WD cooling depends strongly on the critical 
temperature of the nuclear condensate, with higher temperatures corresponding to more relevant neutrino emission[ll]. 
These considerations motivate the detailed study of the thermodynamics of such a nuclear condensate. 

For calculational simplicity, we will work in the regime of stiff electrons, so that they simply screen the Coulomb 
interaction with a screening mass. In this light, our investigation may be seen as an investigation of nonrclativistic 
charged spin-0 bosons interacting via a screened Coulomb (that is, Yukawa) interaction. Surprisingly, we find that 
the system is significantly more complex than expected and that we can merely set an upper bound for the first-order 
transition temperature: Tbec < Tc°\ We conjecture that this low first-order transition temperature is accompanied 
by an unforeseen second-order transition at Tc°^ . 

This paper is organized as follows: in Section II we discuss this model in detail and calculate its one-loop effective 
potential. In Section III we establish the phase diagram described by this model and investigate its properties 
analytically and numerically. In Section IV we demonstrate that the condensed phase is globally disfavored anywhere 
the usual uncondensed phase exists, and to resolve this puzzle conjecture that the phase transitions that this system 
undergoes are more complicated than previously appreciated. Finally, in Section V we make some remarks about the 
phenomenological relevance of nuclear condensates and discuss priorities for more deeply understanding this model. 



II. MODEL, QUASIPARTICLES AND EFFECTIVE POTENTIAL 

At the densities we are considering our system can be described by the (Euclidean space) Lagrangian 



i> + -^F llv F ilv + Cgauge + V{ D ^1 H + m + VeJo)r], (1) 



where if> is the spin-0 helium nucleus field with charge Ze — 2e so that D^ip = d^ip — iZeA^ip, r\ represents the 
usual electron, with D a r\ = d u rj + ieA^rj, A u is the photon field, M (m) is the helium nucleus (electron) mass and 
fi (p e ) is the chemical potential for the nuclei (electrons). C gauge is the gauge-fixing action required for perturbative 
calculations. 

The Lagrangian in (1) does not contain the nuclear force between ions. We omit this force because in the regime 
we are considering the Coulomb repulsion prevents two nuclei from approaching one another to distances comparable 
to the nuclear size and, consequently, the nuclear force between them is inoperative: the nuclei interact with each 
other (and with electrons) only through the electromagnetic force. At the densities we consider (n < 10 s g/cm 3 ) the 
nuclei are non-relativistic while the electrons may or may not be relativistic. A chemical potential for the electron is 
included and chosen so that the charge density of electrons equals that of the nuclei, ensuring charge neutrality. 

Despite its apparent simplicity the action above describes a tremendous array of phenomena; this is not surprising 
giving the number of parameters present. We will concentrate on the regime described in the introduction where three 
different small parameters can be identified, namely i) ami — I /do, the ratio between the particle distances and the 
Bohr radius, ii) the fine structure constant a and iii) the mass ratio m/M. In order to make progress we will attempt 
a calculation that captures the leading order effect on these three parameters. Sometimes, however, certain effects 
will be proportional to ratios of these parameters. In these circumstances we will consider their numerical value to 
decide which terms to neglect. For instance, we count m/(aM) m 10~ 2 as a small parameter. 

The first step is to integrate out the electrons. The result is, in general, a complicated non-local action for nuclei 
and photons only. Later, we will use only the part of the action quadratic in the fields. So, at leading order in a we 
have 
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L,pA = "0 T ( D - /i - — j ip + -F^F^ + Cg aU g e + -A^Il^A,, - (eA + [i e )n. (2) 



The quadratic term in A a is the one-loop photon polarization tensor due to the electrons. Higher loop corrections 
are suppressed by powers of ami = l/ao and are small for the dense electron plasmas considered here [24]. 

In the density and temperature regime considered here, the electrons are degenerate and we can use the T = form 
of the polarization tensor 



n -^n 
.a« n ^n + te-% P 2 )n J 
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where n and H 1 - are functions of Po,p ■ The form of the polarization tensor and the fact that it is determined by 
two functions follow from the Ward identity p M n M „ = 0. H- 1 will play no role in what follows but n will be essential 
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for our discussion. It can be written as H(po,p) — m 2 f(pom/kp,p/kF) where / — > 1 at small pom/kp,p/kp [24]. 
The zero momentum value of II describes the static screening of Coulomb forces by the cold electron gas. For certain 
values of the momentum II also has a small imaginary part. In this paper we will neglect both the imaginary part 
and the momentum dependence of II, effectively working with a model where the electrons provide a negative charge 
background canceling the ion charge and screening for the Coulomb force. In reality the electron gas leads to the 
"Friedel oscillations" in the screened Coulomb force that will be considered in a further publication. Our calculations 
are thus applicable to a model of spinless bosons interacting through a screened Coulomb (Yukawa) potential. The 
effect of the full momentum dependence of II and the Friedel oscillations on the thermodynamics will be left for a 
later publication. 

We are interested in the possibility that nuclei condense, that is, that the condensate (ip) = v be non-zero, breaking 
the electromagnetic U(l) symmetry spontaneously. Whether this happens or not can be decided by minimizing the 
effective potential V(v) . For this purpose we compute now the one-loop effective potential by using standard methods 
[25]. First, we split the nuclear field into a classical part v and a fluctuating piece x as V' = v + XR + i>Xi- Then 
we expand the action to quadratic order in the fields A^, x, X^ an d perform the gaussian path integral. The effective 
potential is given by 



-/3f d 3 rV(v) 



J DxD^DA^ e - s «"< 



= (dct S qU ad) 
_ „-5trlnS,2„ a <j 

where j3 — 1/T is the inverse temperature. 

The fields Xi X* an d ^-p iiux m ^ ne quadratic part of the action. The unitary gauge-fixing used in Ref. [10], 
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followed by the limit £ — > oo allows us to decouple the fields at quadratic order and therefore simplify many calculations. 
However, the use of the unitary gauge (and R% gauges in general) is known to be problematic at finite temperature 
[26]. The gauge fixing condition depends explicitly on v and with this gauge we would be computing the effective 
potential (a gauge-dependent quantity) at different values of v in different gauges. For this reason we will instead use 
the Coulomb gauge-fixing, 
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(6) 



followed by the £ — > limit. For the one-loop calculations described in this paper the use of Coulomb gauge is a 
modest calculational complication that preempts more complicated conceptual questions. 
The quadratic part of the action is given, in momentum space by 
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The eigenvalues of the matrix in (7) are po — iE p and p a 
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(8) 



We are assuming that translation symmetry is not spontaneously broken and v is position independent. 
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where m 2 A — 4naZ 2 v 2 /M. Eq. (8) gives the dispersion relation for the quasiparticles in the system (after the 
analytic continuation po — > ipo). At zero temperature and up to higher loop corrections we can use the tree value of 
the chemical potential fj, = in the dispersion relation above. In that case, the result in (8) becomes independent of 
the gauge fixing parameter £ (as it should since it is an observable quantity) and it agrees with the dispersion relation 
obtained in using the unitary gauge [10]. 

The computation of the one-loop part of the potential can then proceed in the usual fashion. We introduce 5 
a spurious offset to p 2 , + E 2 , do formal manipulation, and then eliminate 5. 
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where the sum is over the discrete values po = 2nTj, j = 0, ±1, Had we kept the full H(po,p) instead of 
merely II = m 2 we would have encountered cuts in the complex po plane and the effective potential would have a 
more complicated form. We will ignore the temperature-independent piece when exploring this effective potential in 
pursuit of its corresponding phase diagram. 



III. PHASE DIAGRAM 



The minimization of V in relation to v gives us the actual expectation value of the condensate at any given value 
of the chemical potential [i. We would like, however, to have the ion density n fixed in order to neutralize the charge 
of the electrons. We then have to solve simultaneously the pair of equations 
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The first terms in these equations come from the tree level contribution V^ ' = —/jv 2 to the effective potential. In 
the absence of dynamical electrons and screening effects (m s = 0), (10) were derived in Ref. [27]. 

We notice that the dispersion relation dependence on /i, shown in equation (8), is important in deriving these 
relations, even if \x is set to zero afterwards. The non-relativistic limit of the dispersion relation obtained in Ref. [23], 
for instance, differs from ours at finite values of \i and is, consequently, at odds with Ref. [27]. 

To the extent that higher loop contributions to the effective potential are small (10) determine the variation of the 
condensate v = v(T) with the temperature. However, as we can see from (10b), /z cannot be negative if v ^ 0. We are 
led then to expect that the condensate forms at positive values of /i. However, at positive /t, E 2 can be negative for 
some values of p, and thus the effective potential is complex at positive /i. A complex one-loop effective potential is a 
common occurrence and signals an instability [28] . Frequently, this instability is an artifact of the loop expansion and 
can be cured by a resummation of higher loop contributions, as it occurs, for instance, in the relativistic \<f) A model 
at finite temperature [26]. We are unable at this point to identify the necessary resummations needed in our model. 
But, just like in the models discussed in Ref. [26], the un-resummed one-loop potential already carries important 
information about the thermodynamics of our problem. We will now proceed to extract as much information from 
the one-loop effective potential as possible. 

The analysis of (10) is simple in the v = case and found in textbooks. In this case equation (10a) becomes 

^ 1 (11 ) 
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FIG. 1: The red curve is the solution to (13a); it is shown as dashed where the condition /i <C p 2 /2M is violated and 
should not be trusted. The result of the phonon (dotted blue) and plasmon (dotted orange) approximations are also 
shown. The purple dashed line separates the phonon-dominated from the plasmon-dominated regions. 



Following the usual analysis, /i is negative for 

T > T(°) = *L (-JLJ) = ( -^—) — (12) 
M\((3/2)J \2({3/2) 2 J MZ 2 ' 1 ' 

and vanishes at T = Tj° . For T < T'j -' it is impossible to satisfy the equations, which signals the need for a 
non-vanishing condensate. The line v(T) = is shown in Fig. 1. 

In order to consider non-zero values of v our strategy will be to neglect \i as compared to p 2 /2M in (10) obtaining 
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with E 2 = (p 2 /2M) 2 + p 2 m 2 A / \p 2 + m 2 ), which obviates the complex effective potential problem. We then solve 
(13) and carefully verify the regime of validity of the /i <C p 2 /2M approximation. We will solve the equations both 
numerically and, in some limits, analytically. The numerical solution of (13a) and ao/l = 35 (corresponding to a 
density of p = 4.6 x 10 5 g/cm 3 is shown as the thick line in Fig. 1. 

Fig. 1 has the typical shape of a first order phase transition. In the usual scenario, for values of temperatures 
where three solutions exist (T^ < T < ST^ in Fig. 1 ), two are locally stable and the middle one is unstable. The 

v = solution is stable for temperatures higher than a critical value T c > rj ^ and the v ^ is only metastable. 
At temperatures lower than T c the roles between the v = and v ^ are reversed. Thus, the condensate v jumps 
discontinuously to zero as the temperature is increased past T = T c . As we will see in the next section, the curve 
shown in figure Fig. 1 cannot be trusted for all values of v and the situation in our model is more complicated. 

We now assess the validity of the approximation fi p 2 /2M leading to our values of v(T). For any given value of v 
and T we can estimate \i by using (13b). This estimate will be accurate if fj, is indeed negligible compared to p 2 /2M 
but not otherwise. In this sense, the use of (13b) conservatively estimates the range of validity of the fj, <C p 2 /2M 
approximation. The estimate of p 2 /2M is a little trickier. Ordinarily, the value of p 2 /2M could be estimated from 
the knowledge of the typical value of p contributing to the integral in (13a). The integrand in (13a), however, has a 
double hump structure dominated by two widely separate scales as shown in Fig. 2. Depending on the values of v, T, 
one or the other hump will dominate the value of the integral. Fortunately, analytical approximations are available 
in these two cases. 

For the lower momentum hump we are in the "phonon region" where the approximations 
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FIG. 2: Integrand in the second equation in (13a) as a function of momentum (full line). The dotted lines 

correspond to the phonon and plasmon approximations. 

are adequate. In this region the dispersion relation approaches that of a phonon [10] 

E p » ^p (15) 

m s 

and the integrals in (13) can be analytically calculated: 

, Mm* „ 

n = v + — — -T , (16a) 
lZniA 

15 M m A 

Equation (16a) determines the condensate v(T) as a function of T. Its value is plotted as the (blue) dotted line in 
Fig. 1. 

For the phonon approximation to (13) to be legitimate it is necessary that the equations (14) be satisfied. The 
integrals in (13) are cutoff by the Boltzman factor e - m Ap/(Tm s ) _ typical value of the momentum is p « 

m s T/mA- Using this value of p, both conditions in (14) become 

T < m A . (17) 

Since v 2 < n, the second condition in (14) follows from the first one. The region excluded by this condition is shown 
as the darker blue area in Fig. 3. We can now compare fi to p 2 /2M. We find 

2 Z 2 tt 3 T 2 
(i « — — — am a — « 1. (18) 
2A1 15 m A 

The region where equation (18) fails is shown in blue in Fig. 3. 

This condition T <C ttia is always satisfied when (18) holds and for the relevant values of the density parameter Ija^. 
We conclude then that, as long as the phonon region dominates the integrals in (13a), we are justified in neglecting 

At higher momentum, on the second hump of the integrand, the "plasmon approximation" is adequate: 

m 2 s < p 2 (19a) 



p 2 Mm 2 , 
2M p 2 



« (19b) 
In this region, the dispersion relation becomes that of a massive excitation (plasmon) 



WQy 2+mi ' (2o) 



the integrals are cutoff by the exponential statistical factor e -^iv 2 /?M)z+m\/T _ e -m A /T e -p 4 /8M 2 mA T ^ Also it turng 
out, (19) implies => T <C \mA and hence we can drop 1 compared to e^ Ep in the integrand. The typical momentum 
is given by 

p 2 « y/8m A TM. (21) 
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FIG. 3: Regions where the phonon (plasmon) approximation is not valid are shown in blue (yellow). The purple 
dashed line separates the phonon-dominated from the plasmon-dominated regions. The red curve is the solution to 
(13a); it is shown as dashed where the condition /i -C p 2 /2M is violated and should not be trusted. 



With the approximations in (19), (13) reduces to 
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We can now verify for which values of v and T the plasmon approximation is valid: 
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The first condition (23a) excludes a tiny region of the w — T plane near T sa where our calculations which neglect the 
temperature independent one-loop effective potential are not valid. The second and third condition are actually very 
similar since the exponential factor e ~ m A/T | g verv small when as T < mj. The areas excluded by these conditions 
are shown in yellow in Fig. 3. 

Finally, the phonon and plasmon region contributions to the integral in (13a) should be compared by taking the 
ratio of (ICa) and (22a). We find 

(n -v 2 )\ , 2 1 / 4 tt 2 m T 7 / 4 

_^ > I phonon * » '»s-t m A /T lnA\ 

("-« 2 )l plasro o n ~12r(5/4)VMmf ' 1 j 

The contour separating the phonon-dominated from the plasmon-dominated regions is shown as a dashed line Fig. 1 
and Fig. 3. For v ~ y/n, the only part of the v = v(T) solution where those approximations make sense, the phonon 
contribution dominates at small temperatures while the plasmon contribution dominates at higher temperatures. 
We have now our full verification of the // = approximation. On the upper branch of the v — v(T) solution up 

to T < 5Tj ' the phonon and the fi ~ approximations are valid. At higher T until T w 7.5Tc°\ the plasmon and 
the fj, « approximations are valid. Beyond that and for the lower branch of v — v(T) < 0.7-^/n our approximations, 
including the neglect of fx, are no longer valid. For this reason we draw the v = v(T) solution as a dashed line in 
Fig. 1 and Fig. 3. 

The numerical examples presented in the figures correspond to a fixed value of the density parameter l/ao = 1/35. 
It turns out that the dependence of v/y/n and other quantities on the density is very mild, particularly if we remember 
that only a relatively narrow range of densities around 10 5 g/cm 3 is phenomcnologically relevant. In fact, in addition 
to the dependence of m s ~ k X J 2 ~ n 1 ^ 6 , only the overall normalization of the integral on (13a) is dependent on n. As 
a rule, however, there is a reduction of the phonon dominated region in the phase diagram as the density is raised, as 
can be seen on (24). 



IV. GLOBAL STABILITY AND A CONJECTURE 



Having identified the locally stable states of the model we now study their global stability. In a certain range of 
temperatures two states, the trivial v = and the condensed u ^ one satisfy the conditions for the local minimization 
of the effective potential. We want now to decide which one is the global minimum. There two states, however, are 
obtained from solving (13a) and, consequently, have the same particle density n. In order to decide which one is the 
stable state we have to compare their free energies defined by 

F(n,T) = V eff (v,T)+ f m, (25) 

where /i and v are the solutions to (13a) and (13b) for given values of T and n. The evaluation of F is particularly 
simple right at T = T c (0) . 

Let us first compute this value for the trivial v = solution. Right at T — Tj? the uncondensed solution has H = 0, 
and we find 
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This is to be compared to the free energy of the v ^ states. The T = Tc , u ^ solution is well within the phonon 
dominated region where the free energy can be computed, with the results in (16), to be 
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F phonon (n,Ti°)) = - C01/3 , B/3 , 0/OWB/2 777F 1 - , 4/3 , Q/0WS/fi A /T7 (28) 
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The free energy of the v — solution is more negative than the free energy of the condensed state by a factor 
proportional to the large parameter (M/m) 3 / 2 . This shows that at all temperatures where the v = state exists it 
is it, and not the condensed v ^ state, which is the globally stable state of the system. A numerical calculation of 
the free energy using the solution of (13a) and (13b) confirms this result and shows that the difference in free energy 

between the two competing states increases as the temperature is further increased past T = Tj ^ . 

This result is at odds with the standard picture of a first order transition. In the usual case, one has -Fphonon < F v= q 

for all T <T C where the critical temperature T c is greater than J 1 ] -* . Even though we were unable to compute F(v, n, T) 
for arbitrary v, it is easy to see that no function F(v,n,T) could have a set of local and global minima as shown in 

Fig. 1 while also globally favoring the v = state all the way down to . To understand this impossibility, consider 
the shape of this purported function around T ~ T'j '' as T is increased. A new minimum at v — is supposed to 
appear in addition to the non-trivial one with v ^ and immediately become the global minimum of the function. 
This is clearly impossible. 

The only way out of this inconsistency is to assume that the state with small v beats the superconducting state 
uM'/iiat temperatures smaller than T c . This can occur if the actual curve v = v(T) has the shape shown in the left 
panel of Fig. 4. In that case a first order transition occurs at a lower smaller than tJ - 1 followed by a second order 
transition. The shape of the free energy as a function of V at different temperatures indicated on the right panel is 
sketched on the right panel of Fig. 4. The extra knee in the v = v(T) curve can occur in the region of the v — T plane 
where our calculation is not under control. The part of our calculation that is under good theoretical control forces 
us to believe in this more exotic possibility 
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FIG. 4: The left panel shows a schematic of our conjectured phase diagram for the condensate v as a function of T. 
The red curve in the left panel shows the minima in F(v, n, T) for a fixed n. The green curves in the right panel are 
cross-sections of F(v, n, T) at fixed T, with the sub-figures in the right panel corresponding to the slices in the left 
panel with the same label. At temperature a only the condensed solution exists. At temperature b the free energy 
develops an inflection point which turns into one minima and one maxima. The first-order phase transition happens 
at temperature c when the wells become equally deep. By temperature d (ie. T^ -*), the well with small v is already 
favored, and there is a second-order transition to v — 0. At higher temperatures the minimum vanishes, shown 

in (e) and (f). 



V. CONCLUSION 



Despite the difficulty posed by the fact that the effective potential is complex for positive chemical potential we 
were able to extract some physical consequences from our one-loop calculation. The main one is the indication of 
a sequence of a first order transition followed by a second order phase transition. This conclusion was developed 
by looking at the free energy values computed within the range of validity of the approximations performed and, 
as such, is quite robust. On the other hand we were not able to compute the critical temperature where the stable 
(superconducting) and metastable (normal) states trade places. But if the double transition conjecture described 
above is correct, this temperature is below the free boson critical temperature Tc°\ 

We also found that the superconducting phase exists, as a metastable state, for temperatures up to about 8 times 
Tc°\ This conclusion agrees with that in Ref. [23]; this is not entirely a coincidence. Contrary to the present paper, 
Ref. [23] analyzes the unscreened model (m s = 0). But the High T region of the v = v(T) curve is in the plasmon- 
dominated region where the screening is not important. In addition, the methodological differences between this 
paper and in Ref. [23] lead to a change in (10a) that is numerically small and, just like here, the [i dependence of the 
dispersion relation is neglected, albeit for different reasons. 

Although the model we analyzed (bosons interacting through a screened Coulomb (Yukawa) potential) is interesting 
on its own merits, applications to high density physics require a proper treatment of the effects of a dynamical electron 
background. Technically the main effect is the inclusion of the contribution of the cuts of II(po,p) in the effective 
potential calculation. Physically they correspond to the fact that the actual force between two bosons presents an 
oscillating component (Friedel oscillations) [29-31]. A detailed of the influence of these effects on the thermodynamics 
of the system will be left for a future publication. 

We have also not discussed the metastability of the superconducting state in a quantitative fashion. In particular, 
we have not estimated its lifetime. This is due to the fact that a proper estimate would us require to compute 
the free energy F as a function v in order to understand the size of the potential barrier separating normal from 
superconducting phases. Until a deeper understanding of the resummations needed to make sense of the one-loop 
results is achieved, this calculation is impossible. In fact, all the results discussed here as a well as a confirmation 
or falsification of our conjectured phase diagram hinge on an understanding of this resummation and that should be 
viewed as the number one priority for further progress in this topic. 

Finally, we can use the results of this paper to assess where the idea of nuclear condensates in dense matter stands. 
Recall that, for the existence of a intermediate temperature regime where the nuclear condensate can exist it is 
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necessary that the crystallization temperature be smaller than the condensation temperature. If one is willing to 
consider metastable states, our estimate that the superconducting state extends up to 8 T'J '' is similar to the 
hypothesis made in [5]. A such, the estimate that, at densities around 10 5 <?/cm 3 relevant for white dwarf physics a 
nuclear condensate should exist, stands unaltered. Only an estimate of the decay time of the false superconducting 
ground state can decide whether the inclusion of the metastable state is appropriate but, considering the extreme 
slow evolution of white dwarfs and the fact that they start out at high temperatures, suggest the the metastable state 
is irrelevant. In that case, only at much higher densities (p > 2.4 x 10 7 g/cm 3 ) and temperatures (T > 10 6 K) can 
the nuclear condensate exist. 
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